% 给定飞机参数
initial_x_plane = 0; % 米
initial_y_plane = 0; % 米
initial_vx_plane = 100; % 米/秒
initial_vy_plane = 0; % 米/秒
gravity = 9.81; % 米/秒^2
time_plane = 0:0.1:100; % 秒

% 给定轮船参数
initial_x_ship = 0; % 米
initial_y_ship = 0; % 米
initial_vx_ship = 10; % 米/秒
initial_vy_ship = 0; % 米/秒
time_ship = 0:0.1:100; % 秒

% 计算飞机的运动
x_plane = initial_x_plane + initial_vx_plane * time_plane;
y_plane = initial_y_plane + initial_vy_plane * time_plane - 0.5 * gravity * time_plane.^2;
v_plane = sqrt(initial_vx_plane^2 + (initial_vy_plane - gravity * time_plane).^2);

% 计算轮船的运动
x_ship = initial_x_ship + initial_vx_ship * time_ship;
y_ship = initial_y_ship + initial_vy_ship * time_ship;
v_ship = sqrt(initial_vx_ship^2 + initial_vy_ship^2);

% 绘制曲线图
figure;
subplot(3, 1, 1);
plot(time_plane, x_plane, 'b-', time_ship, x_ship, 'r--');
xlabel('时间 (秒)');
ylabel('水平位置 (米)');
title('水平位置-时间曲线');
legend('飞机', '轮船');

subplot(3, 1, 2);
plot(time_plane, y_plane, 'b-', time_ship, y_ship, 'r--');
xlabel('时间 (秒)');
ylabel('垂直位置 (米)');
title('垂直位置-时间曲线');
legend('飞机', '轮船');

subplot(3, 1, 3);
plot(time_plane, v_plane, 'b-', time_ship, v_ship, 'r--');
xlabel('时间 (秒)');
ylabel('速度 (米/秒)');
title('速度-时间曲线');
legend('飞机', '轮船');